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Abstract. We present Quantum Monte-Gar lo simulations of an exchange-anisotropic spin- 1/2 
Heisenberg model on a square lattice with nearest and next-nearest neighbor interactions. The 
ground state phase diagram shows two classical magnetically ordered phases for dominating 
antiferromagnetic S'^ -interactions and for large quantum fluctuations a ferromagnetic order in 
the x-y plane. In between a finite region is detected where neither classical nor quantum 
mechanical order, e.g. long-ranged dimer correlations, are found. 



1. Introduction 

The investigation of frustrated spin models has become a rather active field in the past decades 
[IHS] . In particular the appearance of quantum critical behavior and emergence of new quantum 
phases in these systems has connected this field of low temperature physics directly with the 
research on high-Tc superconductors and topological insulators \^r&, and references therein]. 
The frustration induced suppression of classical order especially for systems in two dimensions 
and the corresponding enhancement of quantum fluctuations yield an interesting competition 
between ordered ground states ~ like staggered dimer phases - and quantum disordered states 
as the spin liquid state [SHIT]. 

The analysis of such frustrated models with Quantum Monte-Carlo (QMC) simulations is 
often complicated due to the sign problem which limits the ability of the algorithm for competing 
quantum fluctuations. This is one reason why in the present work we choose the quantum 
fluctuations to be non- frustrating (i.e., the S'^^-interactions are ferromagnetic). In addition 
such a spin model can be mapped exactly onto a model of hard-core bosons with repulsive 
interactions [12^113). These bosonic models can be analyzed in optical lattices [14j and also 
frustrated systems have been realized recently [15j . 

In the following section we will introduce the model and describe the expected ground states 
before explaining our QMC results in Sec. [3l We conclude with a short discussion in Sec. HI 



2. Model 

The model is given by summing over all interactions on nearest neighbor (NN) and next-nearest 
neighbor (NNN) bonds between spin-1/2 operators on a L x L square lattice: 



H = Y^hiStSj + H.c.) + Vi S^S] + UStS; + H.c.) + V2 S^S] 

NN NNN 



(1) 
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Figure 2. Left: The spin stiffness inside 
the disordered region at = —0.1 Vk and 
V2 = 0.51 Vi only shows a small signal for 
an intermediate temperature on the 12 x 12 
lattice. For a larger lattice no finite signal is 
found. Right: The energies are saturated to 
their ground state value for T < 0.05 Vi. 
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Figure 1. Ground state phase diagram 
of the anisotropic Heisenberg model on the 
frustrated square lattice (data from our 
previous work [Hj). For small fluctuations \tk\ 
the classical antiferromagnetic phases survive 
and for large values a ferromagnetic order 
in the x-y plane is stabilized. For the 
intermediate region we find no finite signal of 
any order parameter. 

The tfc < couple the 5^^-direction ferromagnetically and the > are chosen 
antiferromagnetic and, hence, introduce frustration in the model. The nearest neighbor S^- 
coupling favors a Neel state which has the lowest energy in the classical Ising case (t^ = 0) 
for values V2 < Vi/2 and competes with a collinear configuration, i.e., parallel spins in one 
direction and antiparallel spins in the other direction, at the critical point V2 = Vi/2. Classically 
the ground state is highly degenerate at this point [16^ 117]. For large quantum fluctuations 
fa 1/2) a ferromagnetic alignment of all spins in the x-y plane is expected. This corresponds 
to a superfluid state of the above mentioned hard-core bosons p^fTO]. 

In the vicinity of the critical point V2 = ^i/2 and for non-zero fluctuations the emergence 
of quantum ordered phases is possible since quantum fluctuations play a crucial role in low 
dimensional systems. This was discussed for similar models, e.g. in [9]. These phases consist 
of quantum mechanically entangled states coupling two or more spins on neighboring sites. 
Examples are staggered or columnar configurations of dimers on nearest-neighbor bonds. 
However, for the present model these dimers are not given by the spin singlet state of two 
spins but due to the exchange anisotropy by the = triplet state. 



3. Results 

The phase diagram of the model ([1]) is given in Fig. [1] in terms of the strength of fluctuations 
—ti/Vi = —t2/V2 and frustration V2/V1. The phase transition lines were derived in Ref. [11] 
from finite-temperature QMC simulations using an implementation of the Stochastic Series 
Expansion |20y21j which is part of the ALPS project |22H24] . Furthermore, an exchange Monte- 
Carlo update |25H28j was used to overcome thermalization problems in the vicinity of the critical 
point V2 = Vi/2. 

To detect the different magnetic phases we calculated the structure factor which is given 
by the Fourier transformed spin-spin correlation function for the wave vectors q = {it, vr) (Neel 
order) and q = (vr, 0), (0, vr) (collinear order). The ferromagnetic phase - corresponding to 
superfluid order in an analogue bosonic model - can be detected by calculating the spin stiffness 
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Figure 3. For tk = -0.1 14 and V2 = 0.51 Vi , 

the complete structure factor of a 20 x 20 
lattice is shown and only a vanishingly small 

signal at q = (tt, tt) (Neel type) is identified Figure 4. The correlation of spins and dimers 
- probably due to finite-size effects and the ^jt^ respect to the top left site and adjacent 
finite temperature (T = 0.02 Vi). bond are shown in a color coded lattice for 

tk = -0.12 Vfc and V2 = 0.5 ViatT = 0.01 Vi. 
No sign of long-ranged order is visible neither 
in the spin-spin correlation function nor in any 
configuration of parallel dimers or plaquettes. 

(or superfluid density). In the QMC simulations the estimator for this observable is connected 
with the winding number |29j . 

As expected the classical antiferromagnetic states survive for small \tk\ and a direct transition 
between both states is found. However, above the critical point we find for \tk\ > 0.1 Vk 
an opening of a region without classical order in the phase diagram. We performed careful 
calculations of the spin stiffness in this region and found it to be suppressed for lattice sizes 
L > 16 and low temperatures T < 0.05 Vi (Fig. [2]). Only for larger values of \tk\ > 0.175 Vk the 
ferromagnetic order parameter shows a finite and size-independent signal [TTj. 

For a set of parameters in the shaded region of the phase diagram in Fig. [1] we calculated 
further observables. In Fig. [3] we show as an example the complete structure factor of a 20 x 20 
lattice in a color-coded map. Apart from a small enhancement of the values for q ~ (vr, tt) which 
occurs due to finite-size effects and thermal fiuctuations (T = 0.02 Vi) we find no evidence for 
any magnetic order in the 5^-components of commensurate or incommensurate type. 

Since no classical order is found, the emergence of quantum ordered phases as the columnar 
or staggered configuration of dimers was also tested. Therefore, we implemented an estimator 
for the correlation of four spins residing on two different nearest-neighbor bonds [5D]. The 
result of such a calculation is shown in Fig. H] for a 16 x 16 lattice and parameters inside the 
unidentified region. The correlation of parallel and perpendicular dimers on nearest-neighbor 
bonds is plotted in a color-coded lattice which refiects the relative distance between the top 
left dimer and a second dimer. The strength of the correlation is normalized to the value of 
the self-correlation of the reference bond. We also include the results of spin-spin correlation 
functions on the lattice sites. We cannot identify any long-ranged order in the dimer or spin 
correlations and, hence, conclude that for a small region of the phase diagram the system is in 
a quantum disordered ground state. 



4. Discussion 

We performed QMC simulations for an anisotropic Heisenberg model on a square lattice with 
competing interactions. The phase diagram was derived and a finite region without classical 
order was identified. An additional analysis of higher order correlation functions also ruled 
out the appearance of dimer ordered phases. In our recent work we also presented exact 
diagonalization results for a set of parameters inside the disordered region and found a finite gap 
in the energy spectrum which is compatible with the rapid decay of the spin-spin correlation 
functions shown in Fig. [H We conclude from our simulations that we have found a quantum 
mechanically disordered ground state with a finite gap which can be referred to as a gapped 
spin liquid. 

Investigations of a similar model on the honeycomb lattice with additional third-nearest 
neighbor interactions are in progress. This model is known to exhibit emergent quantum phases 
for the exchange- isotropic case [3TH33] . 
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